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The Kondo temperature of a magnetic impurity in a weakly disordered metal is distributed 
due to the randomness in the local exchange coupling, and the local electronic density of 
states (LDOS). We show that in a closed, phase coherent metal particle the resulting distri- 
but ion of T k is strongly asymmetric and scales with the mean level spacing A. Its width is 
2yJ d^/ (3Tk, where A is the mean level spacing, and /3 = 1, 2, with, without time reversal sym- 
metry, respectively. Increasing the density of magnetic impurities, the distribution of Kondo 
temperatures (DKT) is found to become more narrow. Corrections to these results due to 
Anderson localisation are discussed. 



1 Introduction 

A local magnetic impurity is known to change the ground state of a Fermi liquid due to correla- 
tions created by the exchange interaction between its localised spin and the delocalised electrons. 
As a result, the magnetic impurity spin is screened at zero temperature by the formation of the 
Kondo singlet with the conduction electrons. Disorder is affecting the formation of the Kondo 
singlet in various ways. The distribution of exchange couplings, due to random positioning of 
the magnetic impurity in t he host lattice, directly results in a corresponding distribution of 
the Kondo temperature^^!. However, the distribution of LDOS in disordered metals is related 
nontrivially through an integral equation with the Kondo temperature. In a previous work the 
distribution of the Kondo tem pera ture (DKT) has been related directly to the distribution of 
the LDOS at the Fermi energyl^l^. This could have some justification in open systems!^. 

Here we rather consider a magnetic impurity in a closed, phase coherent metal particle, 
where the energy levels Ei are discrete, a Kondo box^'''. Then, the self consistent equation, 
determining the Kondo temperature of a magnetic impurity at position r is^^ 

1 = J{v)Vol. J tanh (^^) . (1) 

The local density of states /o(e,r) is in terms of the local wave function amplitudes tpi{r) of 
eigenstates | Z > of the disordered Hamiltonian H given by 

p{e,r)=J2\^Pl{r)\'5{e-El), (2) 
I 

Thus, Eq. can be rewritten as 
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where x = A/ J, where A = l/uVoL. v is the average density of states which is = m/2TTfi^ 
per spin channel in a 2D electron system (2DES). N = Ep/A, with the Fermi energy Ep and 
K, = Tk/A. xi = Vol. I ipii^r) p is the local eigenfunction probability, si = Ei/A are the Eigen 
energies in units of A as measured relative to Fermi energy Ep. 



2 Distribution of Tr- 



In a disordered, phase coherent metal particle, the electron levels repel each other, and the level 
spacings have the Wigner-Dyson distribution^. In this regime, the square of the wave function 
amplitudes xi and the Eigenenergies s/, are distributed independently, xi obeys the exponential 
Porter-Thomas distribution, which is in the unitary regime, with broken time reversal invariance 
(GUE) and in the time reversal symmetric regime (GOE), given respectively byl^ 

^buE(a;z) = exp(-xz), Pgo¥.{xi) = ^^A_ exp (^-^^ ■ (4) 

The probability to find the 1-th energy level above the Fermi energy at energy si can be obtained 
from an integration over the measure of random matrix theoryl^. For large / ^ 1 these probabil- 
ities are in good approximation Gaussian distributed around the equally spaced level spectrum 
= 2, ... with widths which scale like \/I and depend on the symmetry of the Hamiltonian. 
For broken time reversal invariance (GUE), one obtains: 

Pp{si) « exp (^-c^ifi^^ , (5) 

where C/3 « 7r/4, 4/7r, for (3 = 1, GOE, [5 = 2, GUE, respectively. The distribution of exchange 
couplings XjPj{x) is a function of the random spatial position of the magnetic impurity in the 
host lattice^^J, and can thus be taken to be independent from the distribution of the conduction 
electron probabilities xi. We will therefore take x to be fixed for the moment. The distribution 
of Kondo temperatures Tk can thus in the random matrix theory regime be expressed as, 



PiTK)= l[dxiP{xi)Y[dsiP{si)5 [1 - Y^^Unh ^ |9tkE — ta^l^ 
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Approximating tanh(si/2«;) = 1 for si > k and tanh(si/2K) = si/2k for s; < k, we can 
perform the integrals over x/. We find that the energy levels with si < k enter the distribution 
only through the number rix of the energy levels which are in the interval < si < Tk, 
I = 1,2, ...,nK ■ Thus, for k > 1, the random distribution of the energy levels enters the 
expression for the DKT only through the distribution function of tik- For — > 00 one can 
evaluate the resulting expression in the two limits nx ^ 1 and tlk = 1- One obtains the DKT 
in good approximation in both limits and interpolates the distribution for finite rix, 
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where Kp = ^Jttc^/2, kq = Dexp(— x)/A. Thus, the Kondo temperature Tk has a very asym- 
metric distribution, and scales with the energy level spacing A. While the average value of the 
Kondo temperature is found to be independent of the symmetry of the Hamiltonian, T^, the 
width is larger with unbroken time reversal symmetry, 

(8) 

It vanishes in the infinite volume limit, A ^ 0, suggesting that in a large metallic sample the 
distributio n of the Kondo temperature is only due to local fluctuations of the exchange couplings 
J/A = 1/x'^. When there is time reversal symmetry the energy level repulsion is weaker and 
the tendency to localisation stronger. As a consequence the probability of a vanishing wave 
function at a position r is enhanced, as well as the probability of large wave function splashes. 
This explains that the DKT is wider than in the unitary regime. 
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Figure 1: The distribution of the Kondo temperature Tk as scaled with the level spacing A for the GOE and 

GUE ensemble. 

3 Dependence of the Distribution of the Kondo Temperature on the Concentra- 
tion of magnetic impurities 

Without an external magnetic field, -6 = 0, the crossover between the orthogonal and the 
unitary regime is in a mesoscopic sample of size L, where the phase coherence length L^p exceeds 
L, governed by the parameter^! 

X. = (9) 

with the diffusion constant D = vpl/d, where d is the dimension of diffusion, and / the elastic 
mean free path. The spin scattering rate l/r^ is renormalised by the Kondo correlations with a 
maximal value limited by the unitary limit of the scattering crossection, yielding in 2 dimensions. 
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where um = 1/-R^ is the concentration of magnetic impurities with R the distance between 
them. Thus, the maximal value of the crossover parameter is given in 2D by. 
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where Nm = nML"^ is the number of magnetic impurities in the sample of size L, and g = kpl ^ 1 
is the dimensionless elastic mean free path in the metal. When there are only few magnetic 
impurities, Nm < g, the crossover paramter Xs < 1 is small, and the metallic particle is in 
the orthogonal regime. Increasing the concentration of magnetic impurities, the parameter Xs 
increases. As found in the previous section this is accompanied by an according decrease of the 
width of the DKT. The spin scattering rate scales for A <^ Tk like {A/Tk)'^'^- Thus, 
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and we conclude that the distribution P{Tk) follows for Nm < g < > /A^ the wider 
orthogonal distribution, crossing over to the unitary distribution in the opposite regime, Nm > 
g < > /A^, where < ... > denotes the average taken over the Nm magnetic impurities. 



The superexchange between magnetic impurities competes with the Kondo screening. When 
Jrkky exceeds the Kondo temperature Tk of single magnetic impurities, the spin of magnetic 
impurities is quenched, forming an arra y of classical spins, whose spin scattering rate is by 
factor 1/3 smaller than for free spins'^. Thus the question arises, if the unitary regime of the 
distribution of Kondo temperatures can be reached before the magnetic impurities are quenched 
by the superexchange. While the average of the superexchange interaction, < Jrkky >= 
vanishes its typical value is of the same order as in a clean sample, \J < Jrkky > ~ 

umi^J^ cos{2kFR)^^ and it has a wide, namely lognormal distribution^^. 

In low dimensions d < 2, electrons are localised even for weak disorder g > 1. Recently, it has 
been argued that in this regime the Kondo renormalisation is stopped by the l oca l level spacing 
Ac = l/i'dS,'^, where Ud is the density of states and ^ is the localisation lengt hd. Accordingly, 
it is expected that the distribution of the Kondo temperature converges to a finite width and 
scales with Ac, when size L and the phase coherence length exceeds the localisation length 
^. In general, for finite g, corrections to the distribution of the local density of states, and the 
occurence of anomalously localised states will result in further modifications from the Kondo 
distribution as obtained from random matrix theory above. 

4 Conclusions 

The distribution of the Kondo temperature depends in a phase coherent metal particle on the 
time reversal symmetry. The distribution is strongly asymmetric and scales with the mean level 
spacing A. Its width decrases with increasing concentration of magnetic impurities. 
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